Disorder- Induced Superfluidity in Hardcore Bosons in Two Dimensions 
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We study the effect of disorder on hardcore bosons in two dimensions at the SU(2) symmetric "Heisenberg 
point". We obtain our results with quantum Monte Carlo simulations using the directed loop algorithm. In 
the absence of disorder, the system has no long-range order at finite temperature due to the enhanced sym- 
metry. However, the introduction of a disordered potential, uniformly distributed from —A to A, induces a 
finite-temperature superfluid phase. In particular the diagonal correlation length ^ decreases but the superfluid 
order-parameter correlation function becomes a power-law. A non-monotonic finite-size behavior is noted and 
explained as arising due to ^. We provide evidence that at long distances the effects of a weak disordered po- 
tential can be mimicked by an effective uniform potential with a root-mean- square value: = A/VS- For 
strong disorder, the system becomes a Bose glass insulator. 
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I. INTRODUCTION 

The effect of disorder on two-dimensional strongly corre- 
lated systems has been an important problem in condensed 
matter physics for decades since disorder is inevitable in real 
systems. Experimentally, ^He in porous media,' Josephson 
junction arrays,^ and thin-film superconductors^ are good ex- 
amples. Bosonic models are appropriate in studying these 
systems, even though the latter two are microscopically 
fermionic, because the elementary excitations are bosonic.^ 
When the bosons interact with each other via a strong short- 
range repulsive potential, a hardcore boson Hamiltonian is a 
good starting point. 

The hardcore boson model has been used to study 
zero-temperature phase transitions in a number of 
contexts : 5. 6,7,8,9, lo, ii , 12, 1 3, 1 4 j.jjg simple transition as a function 
of average fiUing j^'^'''* the effect of adding a nearest-neighbor 
repulsive interaction (competition between superfluidity 
and checker-board solid), the transition induced by 
a staggered potential (also superfluid — > checker-board 
solid), and the destruction of superfluidity caused by 
adding disorder (superfluid — > Bose glass). One reason 
for the interest in incorporating nearest-neighbor repulsion 
is that it allows access to a higher symmetry situation: for 
a particular value of the repulsion, the Hamiltonian has a 
full SU(2) symmetry, an enhancement of the usual U(l) 
symmetry associated with boson number. 

A key property at finite temperature is whether the system 
has a non-zero critical temperature or not. In the clean sys- 
tem, for example, at the SU(2) symmetric point (known as 
the Heisenberg point in analogy with quantum antiferromag- 
nets), the critical temperature is zero by the Mermin- Wagner 
theorem.'^ Away from this special case, the symmetry of the 
system is U(l), and one expects that superfluidity is possible 
at a finite temperature; this has been studied numerically by 
several groups.- '^''"'" - 

In this paper, we focus on the effect of adding disorder to 
a system at the SU(2) symmetric point at a finite temperature. 
Though the mean value of the disordered potential is zero, 
it does, of course, break the SU(2) symmetry of the system. 
Thus, we study if disorder can induce superfluidity. 



The Hamiltonian for hardcore bosons with nearest- 
neighbor repulsion which we study here is 



where is the creation (destruction) operator of a boson 

at site z on a square lattice of size Lx L,ni = blbi is the boson 
number operator. Note that if fii = for all sites, the system 
lies at the Heisenberg point [SU(2) symmetry]. The distribu- 
tion of the disorder potential, Hi, is taken to be uniform in 
the interval [—A, A]. We study this model using a quantum 
Monte Carlo method based on the directed loop algorithm. 
This algorithm reduces the autocorrelation time so that we 
may simulate lattices up to size 128 x 128. 

In the weak disorder case, we find strong evidence that 
superfluidity exists at finite temperature: the disordered 
potential added upon the Heisenberg point induces finite- 
temperature superfluidity. However, the system has an inter- 
esting finite size effect causing non-monotonic behavior. We 
understand this effect as arising due to the presence of a large 
but finite density-density coiTelation length, which is natural 
in the presence of a weak uniform potential. When the dis- 
order is strong, the superfluidity eventually disappears; the fi- 
nite compressibility at this transition shows that the insulating 
phase is a Bose glass. 

Our paper is organized as follows. In Sec. II, we show 
results for the superfluid order-parameter susceptibility and 
winding number susceptibility for different disorder strengths, 
and show that we find finite-temperature superfluidity. Sec. Ill 
discusses the results for different correlation functions and ex- 
plains the reason behind the interesting finite-size behavior 
Results for a uniform potential are presented for comparison 
in Sec. IV. Finally, we summarize our results and discuss the 
nature of the insulating phase in Sec. V. 



2 



II. ORDER-PARAMETER SUSCEPTIBILITY AND 
WINDING NUMBER SUSCEPTIBILITY 

To quantify whether the system is superfluid or not, we 
measure two physical observables in our simulations: the 
order-parameter susceptibility and the winding number sus- 
ceptibility. The superfluid order-parameter susceptibility is 
defined as 



where 



T 
4L2 



l/T 



dT 



l/T 



dT'VGp(r,r,r',r') (2) 



Gp(r,r,r',r') = (6(r,r)6t(r' 

and the winding number susceptibility is 

r „ 



b\r,T)b{r',T')) (3) 



(4) 



where Wx is the spatial winding of the boson world lines 
in the x direction, and similarly for Wy. Here b{r, r) is the 
Heisenberg operator defined as e^^6(r)e~^^, and (• • • ) de- 
notes averaging over thermal configurations and disorder real- 
izations. Note that the order-parameter susceptibility is the in- 
tegral of the two-point off-diagonal correlation function, and 
recall that the winding number susceptibility is proportional 
to the superfluid stiffness. 

In the large size limit (L oo), the scaling forms for the 
susceptibilities are 



Xp 



and 



Xv 



B 



C 



(superfluid) 
(insulator) 



(superfluid) 
Dcxp{—L/Lo) (insulator). 



(5) 



(6) 



In the superfluid state, rj < 0.25 and C > 2 where the lim- 
iting values are obtained right at the Berezinskii-Kosterlitz- 
Thouless transition. 

In all simulations, we fix the temperature at T = 1/8 (cor- 
responding to (3 — 8) and e = /3/M = l/4 where M is the 
the number of Trotter-Suzuki^^ slices in the imaginary time 
direction. The lattice size ranges up to 128 x 128. The total 
number of loop updates for the thermal average is 10000, and 
the minimum number of disorder configurations is 100. 

Figure [1] shows the winding number susceptibility (Xw) 
and order-parameter susceptibility (xp) as a function of sys- 
tem size for various disorder strengths. In the clean case, 
A = (Heisenberg point), the winding number suscepti- 
bility decreases monotonically as the system size increases. 
Since there is no ordered phase at finite temperature by the 
Mermin- Wagner theorem and the correlation length is very 
large finite size effects are large. Xp shows an approxi- 
mately power-law form; the fit to Eq. (4) is given in Table 
U For Xw, we fit to the form E — F\ogL expected from 
renormalization group arguments;^" the fit is very good, with 
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FIG. 1 : (color online) Winding number susceptibility (Xw, top, linear 
scale) and order-parameter susceptibility (xp, bottom, log-log scale) 
as a function of system size L for different values of the disorder 
strength A. At both A = and A = 3, the behavior signals the 
absence of superfluidity. However, at the intermediate values A = 
1.0 and 0.1, once the size of the system reaches a certain length scale, 
both the saturation of Xw at a value greater than 2 and the power-law 
increase of Xp with exponent greater than 1.75 signals the presence 
of superfluidity. 



TABLE I: Fitting results for the order-parameter susceptibility (Xp) 
and winding number susceptibility (xw) for A = 0.0, 0.1, and 1.0. 
The fitting formulas are Eqs. Jsj and ([6}; the x^ given is per degree 
of freedom from the fit. 



A 


Xp 

A V 


Xw 

c x' 


included L 


0.0 


0.613(5) 0.218(2) 2.9 




16-128 


0.1 


^0.55 ^0.15 


«3.34 


96,128 


1.0 


0.567(7) 0.151(3) 0.9 


3.572(6) 3.2 


16-128 



E = 4.37(5), F = 0.38(1), and = 1-2 over the range 
Le[16,128]. 

For strong disorder, A — 3.0, as the system size increases, 
the data show that the winding susceptibility decreases expo- 
nentially to zero and the order-parameter susceptibility sat- 
urates. Thus, as expected, superfluidity disappears in the 
strongly disordered case. 
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The most interesting data is for the intermediate disordered 
cases, A = 0.1 and 1.0. For A = 0.1, the winding number sus- 
ceptibility first decreases for small lattices but then increases, 
saturating at a large value of nearly 3.4 for sizes larger than 
L = 96. Similarly, the behavior of the order-parameter sus- 
ceptibiUty seems to change as a function of L. Since we were 
unable to smoothly fit the data with Eqs. ^ and (|6]l, we give 
a crude value, r/ w . 1 5, and C ~ 3 . 34 in Table I obtained from 
using only two points, L — 96 and L = 128, by which size both 
the winding number and order-parameter susceptibilities have 
become stable. For A = 1.0, there is some similar variation 
for the smallest lattices, but the behavior stabilizes for L > 16. 
Reasonable fits yield t] = 0.151(3) and a saturation value of 
Xw of about 3.6. From the values of rj and the saturation of 
Xiv, both of our intermediate disordered cases are superfluid 
at this temperature. 

In summary, we find that there exists a superfluid state when 
A = 0.1 or 1.0, but that superfluidity is absent at A = 3.0 as 
well as in the clean case. Thus, weak disorder induces su- 
perfluidity which is subsequendy destroyed if the disorder be- 
comes too strong. 



III. CORRELATION FUNCTIONS 

The strong and non-monotonic finite size effect in the data 
of Fig.[T]remains to be explained. In this Section, we present 
results for two correlation functions to address this issue. The 
first is the superfluid order-parameter correlation function de- 
fined by 



Cp{r) 



T 
l2 



/ dr dr' VGp(r,r,r',r') 
"'0 J a 



(7) 



The second is the density-density correlation function 



) = J2 dr' J2 Gd(r, r, r', r') (8) 



where 



Grf(r,r,r',r') = (<T(r,r')n(r,T)n(r',T')) 



(9) 



with the definition n(r, t) = 6^'(r, T)b{r, r) — 1/2 and a{r, r') 
is a stagger factor which is 1 if r and r' are on the same sub- 
lattice or —1 otherwise. We take r = x — x', and the sum in 
Eqs. O and (|8]l is constrained such that r is held fixed and 



y = y ■ 



In the clean case, A = (Heisenberg point), we expect 
Cp{r) — 4:Cd{r) due to the SU(2) symmetry of the model. 
The two correlation functions are shown in Fig.|2] Note that 
because of the periodic boundary conditions used in the calcu- 
lation, a correlation function which initially decays must start 
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FIG. 2: Order-parameter correlation function, Cp{r), and density- 
density correlation function, Cd{r), in the clean A = case (log-log 
scale). The error bar is smaller than the size of the symbols, and 
the line is a fit to Eq. Onil. The y-axis is in arbitrary units in order 
to show that the two correlation functions are the same because of 
symmetry. We find the exponent rj = 0.218(2) from Cp(r). 
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FIG. 3: (color online) Order-parameter correlation function [Cp(r), 
upper panel] and density-density correlation function [Cd{r), lower 
panel] for the disordered case A = 0.1. Cp{r) has a power-law form 
for all sizes. In contrast, the behavior of Cd(r) changes from power- 
law decay to exponential decay for L > 32. Thus the diagonal or- 
der disappears for larger systems, leaving the system superfluid. We 
show the error bars of the endpoints. 
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increasing for r > L/2. Thus, here and throughout this paper, 
we show the correlation functions only up to r = L/2. 

It is well known that when A = the system develops a 
finite but large correlation length 3.1 any small temperature. 
When L ^ ^0 it is difficult to detect this correlation length 
and a power-law behavior should be a good description of the 
data. Hence we fit our data for Cp{r) to the form A[l/r^' + 
1/{L - r)**] and find that 77 = 0.218(2), which is consistent 
with the fitting of Xp in Table I] The small wiggle in Cd{r) 
arises from the fact that our Monte Carlo update is optimized 
for measuring the order-parameter susceptibility so that the 
error in Cd is larger than that in Cp. 

Now we turn our attention to the disordered case and focus 
on A = 0.1 for which there is a peculiar dip in the winding 
number susceptibility in Fig.[T| First note that in the presence 
of disorder, the symmetry of the system is affected so that the 
two correlation functions are no longer identical. In Fig.|3] we 
show both the order-parameter correlation function and the 
density-density correlation function for L = 24, 32, 48, and 
64. Thirty disorder realizations are used for each lattice size. 

Fig. 13 shows a change in behavior at L ^ 30. For smaller 
systems, both correlation functions appear power-law. But 
for larger systems, while Cp(r) is truly power-law, Cd{r) be- 
comes exponential - a gap develops in the density-density 
channel. By fitting Cd{r) to the form 



(10) 



we extract ^. We find ^i^s = 17(2) and ^1^4 = 19(1). 
These two values suggest that in the thermodynamic limit ^ 
is roughly 20. 

One can understand the dip in the winding number sus- 
ceptibility and the bump in the order-parameter susceptibil- 
ity (Fig.[TJ with the following scenario: When the system is 
smaller than ^, the system thinks the full symmetry is still 
present, and so the winding number decreases as for A = 0. 
But as the system size grows larger than ^, the density-density 
order is no longer long range, leaving the system in a super- 
fluid state in the thermodynamic limit. As one increases the 
strength of disorder, the length scale ^ becomes small so that 
even a small system shows superfluidity, as for A = 1 .0. Even- 
tually disorder not only destroys the density-density long- 
range order but also destroys the superfluid long-range order 
(A = 3.0). So we conclude that the weak disorder regime has 
severe finite size effects but that disorder enhances superfluid- 
ity before the system finally enters the strong disorder regime. 
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FIG. 4: (color online) For different values of the uniform chemical 
potential /i, the winding number susceptibility, Xw, as a function of 
system size. The behavior is very similar to that in the disordered 
case Fig.[T] no superfluidity at either zero or large fi, but substantial 
superfluidity for intermediate values of p. 



of particles. The finite-temperature phase diagram was stud- 
ied in the anisotropic spin- 1/2 XXZ model.^'^' Here we study 
the isotropic case for comparison with the disordered results 
above and especially focus on the correlation functions. 

The winding number susceptibility, shown in Fig.lH dis- 
plays behavior very similar to that in the disordered case. The 
difference is that even though Xw decreases when fi = 0.1, the 
value itself is larger than in the clean (A = fj, = 0) case. We 
learn that for small sizes weak disordered potential suppresses 
Xw compared to the clean case while a weak uniform potential 
enhances it. The correlation functions at /i = 0.1 in Fig. |5] are 
also similar to those in Fig. [3] Using Eq. ( fTOb to fit, we find 
^1^32 = 13(1), ^^48 - 13.2(5), and Ci^i = 12.8(5). These 
are all in reasonable accord with each other suggesting that ^ 
is roughly 13 in the thermodynamic limit 

For uniform potentials, the hardcore boson model given in 
Eq. ([T]) is equivalent to a quantum spin-half antifeiTomagnet 
in a uniform magnetic field. At zero magnetic field, it is well 
known that the long wavelength physics of the system is de- 
scribed by a non-linear 0(3) sigma model.^^ In particular an 
exact formula for the correlation length at small temper- 
atures has also been derived using this connection. Using 
these ideas we can write down the effective field theory that 
must describe the physics of our model for small /i as long as 
^oAi ^ 1- The action of this theory is given by 



IV. UNIFORM POTENTIAL 



dyS-dyS- fi^Sl+S^)] (11) 



It is instructive to compare the effect of a disordered poten- 
tial with that of uniform potential, /i^ = /i for all i in Eq. ([T]). 
Since a nonzero uniform potential explicitly breaks the sym- 
metry of the Hamiltonian, it is well known that the system is 
superfluid at T = until the uniform potential reaches a criti- 
cal value, beyond which the system enters an insulating phase 
because the hopping is impeded by the the large occupation 



where S = {Sx, Sy, Sz) is a unit three vector The density- 
density coiTelations are mapped into the correlations in Sz, 
while the order parameter coiTelations are mapped into coiTe- 
lations in Sx ■ Expanding the theory about its minimum then 
shows that the fluctuations in Sz have a correlation length 
^ cx I//1 while correlations in Sx are long range. Indeed we 
find that the density-density correlation lengths ^ = 12.8(5), 
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FIG. 5: (color online) For the clean case with uniform potential /i = 
0.1, the order-parameter correlation function (Cp(r), upper panel) 
and density-density correlation function (Cd(r), lower panel). The 
behavior of Cp(r) is power law, while Cd{r) becomes exponential 
when it reaches L = 48. We show the error bars for the endpoints. 



5.8(5), 4.5(5) at /.t = 0.1, 0.2, and 0.3, respectively (for 
L — 64), agree with this prediction. 



V. DISCUSSION 

In this paper, we have studied the presence or absence of 
two types of long-range order at finite temperature. When the 
system has SU(2) symmetry, the superfluid order-parameter 
correlation function and the density-density correlation func- 
tion behave identically, and there is no long-range order We 
tune away from the SU(2) symmetric point either with a uni- 
form potential or with a disordered potential which has zero 
mean. In both cases, we show that the density-density corre- 
lation function decays exponentially, while the super-fluid or- 
der parameter correlation function decays as a power-law, im- 
plying that superfluidity exists for sufficiently large systems. 
For sufficiently large potentials of either type, all correlations 
again become short range and the superfluidity disappears. 

We can easily understand the existence of superfluidity for 
small uniform potentials using an effective field theory; it 




FIG. 6: Compressibility upon applying either a uniform chemical 
potential (left panel) or disordered potential with zero mean (right 
panel). The size of the system is 64 x 64, and the number of disorder 
realizations is 10 in the disordered case. Note that in the disordered 
case we expect 1 < Ac < 3. Vanishing compressibility indicates that 
the non- superfluid phase is an incompressible Mott insulator, while 
finite compressibility indicates a "Bose glass" phase. 



predicts that the density -density correlation length follows 
^ ~ If we assume that the disordered potential can be 
replaced by an effective uniform potential equal to the root- 
mean-square value, we get //off = A/\/3. Using our results 
for the uniform potential, we then predict roughly ^ ~ 22(1) 
for A = 0.1 which agrees very well with our results in the dis- 
ordered case. Thus, we think that a weak disordered potential 
indeed behaves very much like a uniform chemical potential 
at long distances. 

To understand the phase for large /i or A, we measure the 
compressibility as a function of /i or A. The compressibility 
is defined as 



1 



TL2 



[^n(r,T) - {^^n(r,- 



(12) 



where n(r,T) = {r , T)b{r , t) is the boson number at site 
(r,T). 

In Fig. 6, we see a striking difference between the com- 
pressibility of a clean and disordered system. As the uniform 
chemical potential /i increases, the compressibility becomes 
zero at about iic ^ 4. In contrast, the compressibility in the 
disordered case remains non-zero, even though we go well 
beyond Ac which is between 1 and 3. The uniform potential 
increases the boson density until the system becomes an in- 
compressible Mott insulator with density fixed at 1, since it 
is hardcore model. But in the disordered case, because the 
mean potential is 0, the system remains at half filling, allow- 
ing density fluctuations to persist. Nonetheless, the system is 
no longer superfluid; it enters the "Bose glass" phase."* 

We close by emphasizing that finite size effects in this sys- 
tem are more dramatic when disorder is weak: one must aver- 
age over distances greater than ^ (the density-density correla- 
tion length) before the effects of disorder will break the SU(2) 
symmetry. For smaller system sizes disorder suppresses su- 
perfluidity, as measured through Xw, as compared to the clean 
case. In the attractive Hubbard-type fermionic model we re- 
cently studied,^** with nearest-neighbor repulsive interaction 
of fermionic pairs, we find similar finite-size behavior of the 
winding number susceptibility when disorder is weak.^^ If 
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the low-energy theory of this fermionic model is bosonic, as ence Foundation Grant DMR-0506953. 
seems to be the case, behavior similar to that described here 
may occur. For disordered fermion systems, deeper investiga- 
tion is needed into the true nature of the effects of disorder. 
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